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Abstract 

in 

New algebraic approach to analytical calculations of D-dimensional inte- 
grals for multi-loop Feynman diagrams is proposed. We show that the known 
\J2 ■ analytical methods of evaluation of multi-loop Feynman integrals, such as 

. integration by parts and star-triangle relation methods, can be drastically 

^ I simplified by using this algebraic approach. To demonstrate the advantages 

• of the algebraic method of analytical evaluation of multi-loop Feynman dia- 

grams, we calculate ladder diagrams for the massless (f> 3 theory. Using our 
algebraic approach we show that the problem of evaluation of special classes 
of Feynman diagrams reduces to the calculation of the Green functions for 
specific quantum mechanical problems. In particular, the integrals for ladder 
D , massless diagrams in the (ft 3 scalar field theory are given by the Green function 



for the conformal quantum mechanics. 
To the memory of Sergei Gorishnii 1958 - 1988 

1 Introduction 



It is well known that the evaluation of the multiple integrals associated with the 
Feynman diagrams is the main source of physical data in the perturbative quantum 
field theory. Since the number of diagrams grows enormously in higher orders of per- 
turbation theory, the numerical calculations of the integrals for multi-loop Feynman 
diagrams are not sufficient to obtain results with desirable precision. That is why 
analytical calculations of the Feynman diagrams (integrals) start to be important. 

For last few years considerable progress was achieved in analytical calculations 
of multi-loop Feynman integrals (see e.g. [1] [6] and references therein). It is 
interesting that in many cases analytical results for the Feynman integrals are ex- 
pressed in terms of the multiple zeta values and polylogarithms. Note that the 
multiple zeta values and polylogarithms are very interesting and promising subjects 
for investigations in modern mathematics (see e.g. [7], [8]). 

*e-mail address: isaevap@thsunl.jinr.ru 
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The analytical evaluations of the multi-loop Feynman integrals are usually based 
on such powerful methods as the integration by parts [9] and star-triangle (unique- 
ness) relation (see [10], [11] and references therein) methods. These methods have a 
long history. For example, the star-triangle relations have firstly been considered in 
the framework of the conformal field theories [12]. Then it was noticed [13] that the 
star-triangle relation is a kind of the Yang-Baxter equation (see also [14]). In [13], 
this fact was used to calculate the "fishing- net" Feynman diagram (of a sufficiently 
large order) for the four- dimensional (D = 4) 4 theory (as well as "triangle- net" and 
"honey-comb" diagrams for the 6 (D = 3) and 3 (D = 6) theories, respectively). 

In this paper a new algebraic approach to analytical calculations of the massless 
and dimensionally regularized Feynman integrals, e.g., needed for the renormaliza- 
tion group calculations, is developed. In particular, this method is based on using 
the integration by parts and star-triangle (uniqueness) relation methods. The ad- 
vantage of our approach is that we change the manipulations with integrals by the 
manipulations with the algebraic expressions. This drastically simplifies all calcu- 
lations, as it will be demonstrated by some examples. In particular, we calculate 
the integrals for ladder Feynman diagrams arising in the 3 field theory for scalar 
massless particles. The integrals for these diagrams contribute to many important 
physical quantities and have been extensively used in many applications, e.g. in the 
calculations of the conformal four-point correlators in the N = 4 supersymmetric 
Yang-Mills theory [15], [16]. The remarkable fact which we have observed is that 
the evaluation of the integrals for special classes of Feynman diagrams reduces to 
the calculation of the Green functions for specific (integrable) quantum mechanical 
problems. For example, the integrals for ladder massless diagrams in the <f> 3 scalar 
field theory are given by the expansion over the coupling constant of the Green 
function for the D-dimensional conformal quantum mechanics. 

The paper is organized as follows. In Section 2, we outline the basic concepts 
of our operator approach. In Section 3, we explain how the integrals for Feynman 
diagrams can be represented in the operator form. The analytical calculations of 
the integrals for multi-loop ladder diagrams are presented in Section 4. In Section 
4 we also discuss the relation of multi-loop Feynman integrals to Green functions of 
specific quantum mechanical problems. In Conclusion we discuss possible general- 
izations and prospects. 



2 Operator Formalism 

Consider the D-dimensional Euclidean space R D with the coordinates Xi, where 
% = 1,...,D. We denote (xp) = J2i x iPi, x2a — {J2i Xi%i) a and in general the 
parameter a is a complex number. Let g« = q\ and pi = p\ be operators of coordinate 
and momentum, respectively, 

[q k , pj] = i5 kj . (1) 

Introduce the coherent states \x) = \{xi}) and \k) = \{ki}) which diagonalize the 
operators of the coordinate and momentum 

qi\x)=Xi\x), pi\k) = ki\k), (2) 
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and normalize these states as follows: 

(x\k) = exp ^ k i x o) ' / d ° k \k) (k\ = l = J d D x \x) (x\ . (3) 

We also define the inversion operator R 

tf = l, tf = Rq 2D , (x\R=H\, 

(4) 

RqiR = qi/q 2 , K t = RpiR = q 2 p i -2q i (qp) . 
where - := {xi/x 2 }. 

The well-known formula for the Fourier transformation of the function l/k 2/3 can 
be rewritten with the help of eqs. (2) and (3) in the following form: 

1 1 / r(p') \ 

where f3' = D/2 — f3, T(f3) is the Euler gamma-function and f3' ^ 0, —1, —2, . . .. 
Note that a((3) obeys the functional equation a(/3) a(/3') = (2ir)~ D . Formula (5) can 
be interpreted as a definition of the infinite dimensional matrix representation for 
pr 2fi (the matrix "indices" are the coordinates Xi and y^). In this representation the 
operator q 2a is a diagonal matrix 

(x\q 2a \y) = x 2a 5 D (x-y) . (6) 

We call the function l/(x — y) 2a the propagator with the index a. Consider the 
convolution product of two propagators with the indices a and (5: 



/ 



d D z V(a',f3>) 



(X - z) 2a (Z - yfP (X - y)2(°>+P-D/2) ' 



where V(a,(3) = a(a + (3)/ (a(a) a(f3)) . This is nothing but the group relation 
p-2a p-2/3 _ p-2(a +p ) j s wr itten in the matrix form (5). 

Using the definition of the inversion operator R (4) one can deduce the main 
formula 

Rff>R = tf(P+%)f0tf(fi-§) . ( 8 ) 

Then, the group relation (Rp 2(a+ ^ R) = (Rp 2a R) (Rp 2/3 R) is equivalent to the 
identity 

fafifP = fPtfif* , ( 7 = a + /3) . (9) 

which can be represented in the form of the commutativity condition [H a , H-p] = 
for the operators H a = p 2a q 2a . Thus, H a (for all a) generate the commutative set 
of elements in the algebra of functions of pi. Identity (9), represented in the 
matrix form (5), (6), is the famous star-triangle relation (a + j3 — 7) : 

*>* - (10) 



(x - z) 2a ' z^ (z - y) 2 ? (x) 2 ? (x - y)^' (y) 2a ' 
Consider the dilatation operator H := (i/2)(pq + qp) = —H^ which satisfies: 

Hq t = q t (H+l), (11) 
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and generates the si (2) algebra together with the elements q 2 and p 2 : 

[q 2 ,p 2 ]=4H, [H,q 2 ] = 2q 2 , [H,p 2 } = -2p 2 . (12) 

It is known that the special conformal transformation generators IQ (4), the di- 
latation operator H (11), the elements pi and = cjipj — q^pi generate the IP- 
dimensional conformal algebra so(D,2). 

Below the following relations will be important: 

[q 2 , j3 2(Q+1) ] = 4 (a + 1) (H + a) p 2a , (13) 
[<f (a+1) ,p 2 ] = 4 (a + 1) (H - a) q 2a , (14) 
H q 2a = q 2a (H + 2a) , H p 2a = p 2a (H - 2a) . (15) 

These relations are easily deduced from the Heisenberg algebra (1). We also in- 
troduce the notion of degree for the operators which are homogeneous functions 
V>(pi, q_j) of the generators of the Heisenberg algebra. The additive number deg(^) 
is called the degree of the operator ■?/> if [■?/>, H] = deg(^) ■0. In particular, we find 
deg(j)j) = —deg((jj) = 1 in view of (11). 

Using relations (13) and (15) one can deduce the algebraic identity 



4(2 7 -a-/3)p 2Q g 2 V /3 = 
= ^[?, P 2(a+1) ] q 2 "P W ~ ^ P 2a P 2(m) ] • 

In the matrix form (5), (6) this identity looks like 



(16) 



(D - 2a 2 - a 1 - a 3 ) (x\a 1} a 2 , a 3 \y) = 
= a 1 [{x\at, a 2 , a 3 \y) - x 2 (x\af, a 2 , a 3 \y)]+ (17) 
+a 3 [{x\a u a 2 , afi\y) - y 2 (x\a u a 2 , aj 



where af — ctjil, a — —a[, 7 = — a 2 , (3 = —a 3 and we have introduced the concise 
notation for the "vertex" integral 

/ I n= f dPz 1 , 111, , 

{ x\a u a 2 , a 3 \y } j ^ _ ^ ^ ^ _ ^ a ^ a ^ \ x \ ^ ^ p a ' 3 W ■ 

(18) 

Identity (17) is called the integration by parts (or "triangle") rule [9] and plays 
a very important role in almost all analytical calculations of multi-loop Feynman 
integrals. 

The integration by parts rule (16) is a combination of two more fundamental 
relations. The first one is 

4(a - 7 )p 2a e< = 4rP 2a [<7 2(7+1) - f] ~ ~~TJ ^ ' ( 19 ) 

7+1 a + 1 

and the second one is obtained from (19) by the Hermitian conjugation (for real a 
and 7). These relations can also be deduced from (13) - (15). 

Consider the sequence of products of the operators p 2 " 2 *- 1 , q 2oi2k 

i> & : = v(- • • , ...) = • • -P 2ai q 2a2 P 2ai f a4 ■■■ ■ (20) 
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Then the finite difference equation (19) and its conjugated version can be written 
in the form of equations Liip(ai) = where the finite difference operators 

Li = 4(a i - a i+ i) + V; + V m 

generate the algebra with the commutation relations: [Li, L i+k ] = for (k > 3) and 
[Li_! - L u L l+l - L l+2 ] = , [L u L i+2 ] = e (*+i+*+*) L i+1 . 

This W- type algebra has the central element Z = — \ Y Lj which is equal to 
the degree of the operator (20). 

Another set of equations a,i ipa = ipa for the functions (20) follows from the 
star-triangle identity (9). The action of the operators Cj is 

Oi . . oti, a i+l , a i+2 , a i+3 , . . .) = 
= ^(...a,i - a i+ i + a i+2 , a i+2 , a i+1 , a i+l - a i+2 + a i+3 , . . .) . 

The operators (21) define the group S even <g> S dd which is a direct product of two 
symmetric groups generated by the sets of even {cr 2i } and odd {(r 2i+1 } elements, 
respectively. One can show directly that the whole set of the relations Liip — and 
<7fc — if) is closed and combinations of these relations do not lead to new constraints 
on tp (20). 

It is well known [17] that the dimensional regularization procedure requires the 
rule: / d D x/x 2a = for all a ^ D/2. A more precise statement [18] is 

/ dDx ^D, W a) =* to D 6(a), (aeRe) (22) 

where tip = * s ^ ne area °^ the un ^ hypersphere in R D . 

In the framework of the dimensional regularization scheme we extend the def- 
inition of the integral in the left-hand side of (22) for arbitrary complex numbers 

a = |a|e iarg(Q) as 

/ * D *-j^ = *n D 8(\a\), (23) 

where S(\a\) is the radial delta-function. It is clear that: f(a) ■ S(\a\) = /(0) • 5(\ct\) 
for analytic in a = functions f(a). 

The important consequence of the definition (23) is that now one can introduce 
the notion of the trace for the operators (20) 

Tr(^)= Jd D x(x\^\x). (24) 

It follows from (23) that these traces are proportional to the delta-function: 

Tr(^) =7rn D c ai 5(\(3\) , [c^ = J dx(x\^\x)^ (25) 

where j3 = ^deg^^), Xi = Xi/r is a unit vector in R D , d D x = Qo r ' 1 dr dx 
and c a . is a coefficient function which contains the whole information about the 



V,; 



e d>± 
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correlator (x\ ip a \x) = c ai x~ 2 ( D l 2+ @\ The simple algebraic arguments which lead to 
the equation (25) are the following. Note that e t( - H+ ~^ \x) = |e~*x) and for r = e* 
we obtain 

(x\^\x) = (x\e t{H -^^ a e- tiH+ ^ \x) = e" i(2/3+D) (x\fa\x) 
which is consistent with (25). 

Remark 1. In the case deg(-^a) = 0, the operators ip a (20) (which are cut off 
by the conditions = for i < and i > 2k) are expressed as a product of the 
commutative operators H a and their inverse: 

fa = H ai H~l_^ H ai _ a2+a3 ■ ■ ■ H- 2 \_ a2k i H a2k , (deg(H a ) = 0) . (26) 

In this case the problem of calculation of the trace (24) reduces to the spectral 
problem for the commutative operators H a . 

Remark 2. The star-triangle identity (9) can be generalized to 

E,,., n {p 2a htiM) F^htiM (p 2/3 ^,.,„(p)) = 

(27) 

where ik {.) (a = 1,2,3) are any tensors being A;th-order homogeneous polyno- 
mials in fa or pi (deg(/i il ... ifc (p)) = - deg(/i il ... ifc (g)) = k) . 

Remark 3. The following relation holds: 

1 a{ f\*) = TT^ a MW- P8) 



(q - y) 2 / 3 p 2a ' (q - x) 2a p 2/3 ' 

Indeed, the contraction of both sides of (28) with an arbitrary state (z\ and relation 
(5) give the identity. One deduces from (28) that the vector a(a) (q — x) 2a p~ 2a \x) 
is independent of the parameters a and Xi and, thus, should be canceled by the 
operators pi. 

Remark 4. The dual analog R of the inversion operator (4) can be defined as 
RpiR = pi /p 2 , RqiR = p 2 qi-2pi (p q) , 

(29) 

where R 2 = 1, R 1 = Rp 2D and (k\R = 

Remark 5. Certain algebraic identities considered above (e.g. eqs. (13) and (14)) 
are related to each other by the obvious Z 2 symmetry p 2 <-> q 2 , H <-> —H which is 
the well known automorphism of the sl(2) algebra (12). 
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3 The Diagrams 



The Feynman diagrams which will be considered in this paper are graphs with 
vertices connected by lines labeled by numbers (indices). With each vertex we 
associate the point in the D-dimensional space R D while the lines of the graph 
(with index a) are associated with the propagator 



x 



a 



y 



= l/(x-y) 



2a 




The boldface vertices • denote that the corresponding points are integrated over R D . 
These diagrams are called the Feynman diagrams in the configuration space. 

The figures and operator form of integral expressions for the 3-point, 2-point 
diagrams and the tetrahedron vacuum diagram are 





f(an) ■ (x\p- 2 <q- 2a2 p- 2a 'iq- 2a4 p~ 2 <\y} 
Fig.l 



= /(at) ■ (x\p- 2a 'iq~ 2a2 p~ 2a 3q- 2a4 p- 2 <\x) 
Fig.2 

f(an) ■ J d D x (x\p- 2a 'iq- 2a2 p- 2a 3q- 2a4 p- 2a 'sq- 2ae \x} 
Fig.3 

where the factor /(a,) is equal to (27r) 3D a(ai)a(a 3 )a(a 5 ) . 

The two- loop propagator- type diagram on Fig. 2 (which is called in [19] the 
master two-loop diagram) is obtained from the diagram in Fig. 1 in the case x = y. 
The expression for the vacuum "tetrahedron" diagram in Fig.3 is the result of the 
integration of the expression for the master two-loop diagram with the additional 
propagator l/x 2ae . This expression is represented in the form of the trace (24) 





/(act) ■ Tr(p- 2a 'iq~ 2a2 p- 2a 3q- 2a4 p- 2 <q~ 2a6 ) . 



(30) 



The tetrahedral symmetry [18] of this diagram becomes evident if we impose, in 
addition to the cyclic symmetry of the trace {s\ = 1): 



(31) 



Si : Tr(p~ 2a ^~ 2a2 p~ 2a ^~ 2a4 p~ 2a 5q~ 2a6 ) = 
the identity (the reflection of the diagram in Fig.3 with respect to the vertical line; 
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this reflection is equivalent to the permutation of the vertices x and 0; s| = 1): 

(32) 

_ 0,(02)0.(014) rp / fi-2a' 2 g-2ai fi-2a' 3 fi-2oi5 fi-2a' 4 q-2<X6) 

The last identity is deduced from (28). With the help of the transformations si 
and S2 one can permute every two vertices of the tetrahedron on Fig.3. Thus, the 
elements s±, s 2 generate the permutation group S4 (e.g. one can check that the 
element (S1S2) generates the 4th-order symmetry (S1S2) 4 = 1). 

Following the paper [18] we add to (31) and (32) the symmetry (cf. with (21)) 
deduced from the star-triangle equation (sjj = 1): 

s 3 : Tr{p^ 2a ^q^ 2a2 p^ 2al, iq~ 2a4 p~ 2a ^q~ 2ati ) = 

(33) 

= 2 n ^^ _2 ( a i _a 2+a3)g-2a3^-2a2g-2(a2-Q3+a4)jr)-2a5^-2a6^ 

The symmetries si, s 2 , S3 (31) - (33) can be realized as linear transformations of the 
indices oti and, as it was shown in [20], they generate the 1440-dimensional finite 
group S 6 <E> Z2, where S 6 is a group of permutations of 6 objects and Z 2 is a 2-fold 
cyclic group. The most general analytical result for the master two-loop diagram 
on Fig. 2 has been achieved in [21] and [22]. According to this result, the trace (30) 
for a 4 = a 2 — 1 can be expressed in terms of 3F2 hypergeometric series. 

At the end of this section we stress that to write down more complicated diagrams 
in the operator form, we need to extend the Heisenberg algebra (1) to the multi- 
particle case [q^\ pf 1 } = i 5 kj 5 ab , where a, b = 1, 2, . . . are the numbers of particles. 



4 Applications. 

In this Section we demonstrate how our algebraic methods work using the example of 
the L-loop ladder (L-boxes) diagrams. We start with dimensionally and analytically 
regularized massless integrals 



D L (po,PL+i,p;a,P,l) 



d u p k 



R j p? k (Pk- P ) 



n 



1 



b (Pm+l -Pm) 27 " 



(34) 



which correspond to the diagram {x — po, y — Pl+i, z — p): 






Fig. 4 

The dual to this diagram is the ladder L-loop diagram for the massless 3 theory 
(the loops in Fig. 5 correspond to the boldface vertices in Fig. 4): 
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Po-P Pi-P 




Pl+i - P 



Pl+il 



(Pmk =Pm-Pk) 



PO Pi PL+1 

Fig. 5 

This "momentum space" diagram represents the same integral (34), but in another 
graphical form. Using (3) and (5) we obtain for (34) the following representation: 



( L 1 

D L (x, y, z; a, /3, 7) = M] ——, 

\k=o a yik 



■ \y) ■ (35) 



i)J X 'p^o g2a k { q _ z) 2f1 k £ 

Now we simplify the problem by fixing indices of the propagators as = a, 
Pk = P, Ik = 7, and consider the generating function 

00 

D g (x, y, z; a, (3, 7) := 0(7') ^ g L D L (x, y, z; a, (3, 7) = 

L=0 



(x\ [f 1 ' - — 
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y) 



(36) 



q 2a (q-z) 2 ^ 

where g = g/a(j') is the renormalized coupling constant. The right-hand side of 
(36) is the Green function for the two-center operator 



q 2a (q - zfP ' 

This Green function has the following symmetries under the linear and inverse trans- 
formations of the coordinates (x,y,z): 



D g (x,y,z; a,P,i) = D g (y,x,z; a, (3, 7) , 



D g (x -z,y-z, -z; (3, a, 7) = D g (x, y, z; a, (3, 7) , 
1 1 1 
x y z 
Here 1/x := {xi/x 2 } and 



D gz-w{-'-i-} a ,(3,l) = x 21 V 21 D g (x,y,z] a, /3,7) 



5 = 27' -a-(3 = D-{a + (3 + 2-f) . 



(37) 
(38) 
(39) 

(40) 



The proof of (37) and (38) is evident. The symmetry (39) follows from the chain of 
equalities 

/ \ -1 

9 \ 



(x\R 2 [f^' 



q 2a (q - z) 



2/3 



y) = 



-2/3 



q -2(a+{3)( q _ 1)2/3 



R\y)= (41) 



-27 -27 /i 1 j z 2/3 V 

V { x l [ P -q2«(q-lW) 'y 



where equations (4) and (8) where applied. 

Using identities (38) and (39) one can change the indices of the propagators in 
(34). In the case 5 = (see (40)) the combination of the symmetries (38) and (39) 
gives 

D g (x, y, z; a, (3, 7) = x' 2 " 1 y~ 2 "< G Y:P (u, v) (42) 

where 

G,,Au,v) = (u\(p 2 ~<'-^) \v), (43) 

and 



- -2/3 9 x i z i Vi z i 

9z 9 Z / 7\ Ta i u i ~ 9 9 1 v i ~~ 9 9 • 

0(7 ) z zp x A z z y z z A 

The function Gy t p(x, y) is obviously related to the function (36) (if we put there 
j3 — 0, g — g z ), and for the case g z = const from (37), (39) we obtain the symmetry 

G^(u, v) = G^(v, u) = (u 2 v 2 )^ G 7i27 _^(-, -) , (44) 

(note that here we have made the redefinition: 7' — > 7). 

Thus, the problem of analytical evaluation of the integrals (34) for the ladder 
diagrams with the special choice of the line indices (a + (3 + 27 = D) is reduced to 
the problem of calculation of the Green function for the Hamiltonian 

which describes (for 7' = 1) the propagation of the scalar particle in the external 
field of one source fixed in the origin. Now the most interesting cases are: 

1. ) 7' = /3 — 1, (a — 1) which corresponds to the D-dimensional conformal mechan- 
ics; 

2. ) 7' = 1, (3 — D/2 — 1, (a = 3 — y) which is related to the standard D-dimensional 
problem of the dynamics of a charged particle moving in the field of the pointlike 
source. For D = 4 both the cases are equivalent to each other and the Green func- 
tion (43) defines the generating function for the ladder diagrams in the 3 (D = 4) 
field theory. 

Here we calculate the Green function (43) and, therefore, the integral (34) for 
general value of D and for the choice of the indices — D/2 — 1, a — f3 — \ 
(D-dimensional conformal mechanics) with the help of the operator approach. Our 
method is based on the identity: 



1 SK 1 



p 2 -g/q 2 L=0 \ 4 



y (H - 1 + a) L+1 p 2 q 



(45) 



where we denote [. . \ a L = (d^, [...]) . Identity (45) can be proved immedi- 
ately if one acts on both the sides by 7i = p 2 — g/q 2 and uses relations (13)-(15) 
and 

representation for the rational function of H in (45) 
(H-1) (-1)^+1 



(h — 2a)(h — a) L 1 J = 5l j0 , Wh ^ 0. Taking into account the integral 

i in (45) 

J o °° dtt L e ta d t (e^ H -^) 



(H-l + a) L+l L\ 
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the Green function G\^(u,v) (43) is written in the form 

1 



H 



(P 2 - 9z/q 2 ) 



L=0 



(46) 



where for the functions $l(u,v) we obtain the integral representation 



r 

■(u,v) = - / 
Jo 

= -a(l) 



dtt 1 



dtt 1 



(u 



u 



/ -la 

Note that the first function is: <3> = a(l) (u — i>) 2 ( 1_ t) 
and in view of (44) all $^ possess the symmetry 



l \ (u — e t v) 2 , 



(47) 



(as it also follows from (46)) 



The integral (47) for D = 4 reproduces the results of [23] and [24]. For appli- 
cations it is worth having the general expression (47) for D = 4 — 2e and expand 
$l(u, v) over small e: 



&l(u,v) 



4^M^l! $i ( ^ 2) 



(48) 



Here the dimensionless parameters 21,2:2 are defined by the equations: 21 + 2:2 = 
2{uv)/u 2 and 21 2 2 = v 2 /u 2 . The coefficient functions <&£ m the expansion (48) can 
be represented in the following form: 



ln(2 1 2 2 )K(2L-/) ' 



/=0 /' /)' m=0 

where C™ = m! ^j^yi and we denote (k > m): 

F(k — m 



J2(-) m Cri m ( Zl ,z 2 ;2L + l-f) , (49) 



m! (I— m)\ 

Z m (2i, 2 2 ; fc) = 



y- _Vfl 



2i T~ 2 2 



2i 2q 



(^-^2) {n tr=i(Eo 



\k—r. 
n i) \i=l 



n 



For example, we have 



Z (21,22; k) 



r(fc) 

21-22 



(Li fc (zi) - Li fc (2 2 )) , 



Zi(2i, 2 2 ; fc) = (Li fc _i i i(2i, 1) + Lifc-i^zi, f ) - (21 <-> 2 2 )) , 

and in general the functions Z m (zi, z 2 ; k) are linear combinations of multiple poly- 
logarithms 



Lifc-m,m 1 ,..,m I .(w , U>1, • • • , W r ) = ^ 



rao>m>--->nr>0 ""I ' • ' "v 
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2 

with J2 m i — m an d arguments Wi belong to the set {l,z±,Z2, ^, • • •}• The 
first coefficient function in (49) (related to the integral (47) for D — 4) 



2 



^\^z 2 ) = -J— E ( ] / u f L hi'fe) [LW^O - U 2L _ f (z 2 )\ , 

(50) 

has firstly been evaluated in [23]. The next coefficient is 

^\ Zl , z 2 ) = £ QL-^n-W ( Zo ^' n + " Zl ^' Z2 ' n + ^ ' 

We also present the result for the integral (47) in the case u = v (x = y) and 
arbitrary D < 2(L+ 1) which is needed for evaluation of the propagator-type ladder 
diagrams: 

1 ' r(D-2) ^H(t + r>/2-i)f 

For D = 4 we immediately obtain w 2 $l(w, m) ~ C(2£ — 1). 
Remark 1. Relation (14) can be rewritten in the form 

Jo 



p 2 * p 2 4a(H -a -l) 



4a 



-2a J_ 

/ 5 ~o 



In the matrix representation (5), (6) this identity gives the analytical expression for 
the vertex integral (18) for the special choice of the indices (one index is arbitrary) 



(x 1,1 -a, ly = — — / dte**" 1 ^ '- yy lD ' 51 

yi 2 ' '2 |y/ 4aa(l) Jo ( e ^ _ y )2(#-i) v ; 

It reproduces the results presented in [26], [27] since the indices of the lines in the 
left-hand side of (51) are changed by means of the star-triangle symmetries (21). 

Remark 2. The one- and two-loop ladder integrals (34) for the case D = 4 and 
ttfc = Pk — Ik = 1 (which are given by the functions $^°' ) and $ 2 °^ (50)) have been 
found in [25]. This result was generalized in [23] and [24] where (50) for all L and 
the generating function (46) (for the case D = 4 and ctk = fit = Ik = 1) has been 
calculated. The analytical results for 1,2,3-loop (boxes) integrals (34), general value 
of D and special indices of the lines ak,Pk,lk can be found in [4], [27], [28], [29], 
[30] (see also references therein). 

Remark 3. In the case of the conformal indices: a k + (3k + lk + 7fe-i = D (V7c) the 
inversion method (41) reduces the 4-point functions (34), (35) (represented on Fig. 
4) to the 3-point functions. Moreover, the inversion method (41) gives the possibility 
to obtain the duality transformation for the D-dimensional Green function 

l r \ I LA = / 2 2x( a -f )/i I \ I I\ 

K l p 2 <* + g(q 2 )- 2a + E ]y/ y U> K x 1 p 2a + g + E {q 2 )- 2 « 1 y' 
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which relates (for the fixed parameters a — 1, E — 0) the D-dimensional Green 
function for the Hamiltonian Ti = p 2 + g (q 2 )~ 2 to the propagator of a massive {m 2 = 
gfi 2 ) free particle (the dimension parameter /i can be restored by the renormalization 

P ~^ P/f 1 an d Q ~^ A 4 *?) • 

The inversion method (41) in the momentum representation (with the help of 

the operator R (29)) gives the duality transformation for the Green function 

where g and E are dimensionless parameters. For a — 1,D — 3 we have deduced 
the duality transformation for the Green function of the Coulomb problem. 



5 Conclusion 

In this paper the algebraic approach to analytical calculations of the integrals for 
multi-loop Feynman diagrams is developed. Possible modifications and applications 
of the approach are the following. First of all we plan to extend the method to 
the case of massive propagators. It is rather trivial to generalize the integration by 
parts identities (16), (19) to the massive case. The massive generalization of the 
star-triangle identities (9), (27) is unknown (see, however, Sect. 3.2 in the paper 
[3])- 

The algebraic formulation of the Gegenbauer polynomial technique [31], [22] is 
also possible and based on the identity 



1 1 



1 00 r r 

^7 = £<T 2 2 d D zd D w\z)C:(~zw)(w\ + (w^z) 



r n-2a ^ (52) 



where zw = (wz)/(w 2 z 2 )^ and C"(t) are Gegenbauer polynomials. It is tempting 
to apply (52) for the evaluation of the two-loop master diagram on Fig. 2. 

Finally we stress that our methods resemble the methods used in the theory 
of the quantum integrable systems. For example, the star-triangle identity (10) is 
a kind of the Yang-Baxter equation [13], [14] and can be used for the i?-matrix 
formulation [32] of the integrable noncompact Heisenberg spin chain proposed by 
Lipatov [33] for describing the high-energy scattering of hadrons in QCD. So we hope 
that the algebraic approach developed in this paper will be helpful in investigations 
of the Lipatov's model. 
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